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The squeezed limit of tlie three-point function of cosmological perturbations is a powerful dis- 
criminant of different models of the early Universe. We present a conceptually simple and complete 
framework to relate any primordial bispectrum in this limit to late time observables, such as the 
CMB temperature bispectrum and the scale-dependent halo bias. We employ a scries of convenient 
coordinate transformations to capture the loading non-linear effects of cosmological perturbation 
theory on these observables. This makes crucial use of Fermi Normal Coordinates and their con- 
formal generalization, which we introduce here and discuss in detail. As an example, we apply our 
formalism to standard slow-roll single-field inflation. We show explicitly that Maldacena's results 
for the squeezed limits of the scalar bispectrum [proportional to (n^ — 1) in comoving gauge] and the 
tensor-scalar-scalar bispectrum lead to no deviations from a Gaussian universe, except for projec- 
tion effects. In particular, the primordial contributions to the squeezed CMB bispectrum and scale 
dependent halo bias vanish, and there are no "fossil" correlations between long-wavelength tensor 
perturbations and small-scale perturbations. The contributions to observed correlations are then 
only due to projection effects such as gravitational lensing and redshift perturbations. 



I. INTRODUCTION 

In this paper we present a simple and complete frame- 
work to derive the late-time physical observables pro- 
duced by the squeezed limit of primordial three-point 
functions, in particular those produced during inflation. 
These observables include the squeezed limit of the CMB 
bispectrum and the scale-dependent bias for large-scale 
structure tracers which crucially depends on the bispec- 
trum in the squeezed limit [1-3]. The main obstacle one 
encounters is the fact that, in order to calculate the bis- 
pectrum consistently, one has to work to second order in 
perturbation theory throughout. We show how this ob- 
stacle is overcome when one is interested in the squeezed 
limit. We take advantage of the fact that different sets 
of coordinates are convenient in diflFerent cosmological 
epochs and there is a careful choice that makes the com- 
putation easy and transparent. In the remainder of the 
introduction, we present a convenient way to choose co- 
ordinates by discussing the different steps that are sum- 
marized in Fig. 1. 

A first step consists in computing the primordial corre- 
lators generated during the phase of primordial inflation. 
For this purpose, coordinates are used that are comoving 
with respect to some unperturbed universe (see e.g. [4]). 
A second step consists in describing the local late-time 
physical processes, such as for example the decoupling 
of photons during recombination or the formation of a 
dark matter halo, that lead to the generation of some 
light signal that eventually will reach the Earth. Since 
these processes take place in regions that are smaller in 
size than the curvature of the background metric (horizon 
scale H~^), they are best described in Fermi Normal Co- 
ordinates (FNC [5]): this is the unique frame (up to three 
Euler angles) constructed along a timelike geodesic pass- 



ing through a point P in which the metric is Minkowski 
with corrections going as the spatial distance rp from 
the central geodesic squared. In an unperturbed FRW 
Universe, the corrections thus scale as {Hrp)^, while in 
the presence of a metric perturbation h^i, the corrections 
involve second derivatives of h^v (in particular, a con- 
stant or pure gradient metric perturbation is removed on 
all scales). The last step consists in relating the signals 
from the point of emission, for example the last scatter- 
ing surface, or the position of a halo, to the observables 
measured on Earth, where scientists most naturally use 
FNC centered on the Earth's world line. 

In general, calculating the evolution of the bispectrum 
of perturbations in any coordinate frame including modes 
outside the horizon requires a full second-order relativis- 
tic calculation in cosmological perturbation theory, in- 
cluding gravity, baryon physics, radiation transfer, etc. 
This is in general a very complicated problem. How- 
ever, a drastic simplification takes place when we limit 
ourselves to a particular configuration of the bispectrum 
in which two of the momenta, say fci and k2 are much 
larger than the third, fc;;. This is often referred to as the 
squeezed limit with ^3 = Ul being the long mode and 
fci ~ A;2 ~ ks being the short modes. 

Then, as described above, we can make use of the fact 
that in a fc^/Zcs <C 1 expansion, the leading and next-to- 
leading order gravitational effects of the long-wavelength 
perturbation fc^ on the short modes can be removed by 
transforming to the local FNC frame, as was empha- 
sized in [7-10]. The first non-trivial effect comes from 
the second spatial and time derivatives of the perturba- 
tion, which are suppressed by 0(fc|/A:|) (see Sec. Ill for 
a more precise statement). We will ignore terms of this 
order in this paper, as they correspond to actual physi- 
cal interactions, which will need to be treated in detail. 
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Generation Evolution '"^ Emission °^ Observation 

comoving coordinates ^q. (12) ™^ Eq. (38) FNC(emitter) cosmic rulers [6] FNC(observer) 

FIG. 1: Sequence of coordinates employed in the computation of obscrvables predicted by a primordial bispectrum in the 
squeezed limit. The arrows represent the change from one set of coordinates to the next, indicating when this transformation 
is most conveniently performed and with a reference to the relevant equations. Tem and Tobs denote the conformal times at 
which the observed photons were emitted and observed, respectively. 



However, up to this order, in the local FNC frame we can 
simply use linear theory. The non-Gaussian correlation 
will then be captured by averaging over an ensemble of 
FNC, each constructed to remove a different long mode. 

In going from comoving coordinates to FNC, one sub- 
tlety arises from the fact that FNC arc affected not only 
by perturbations but also by the expansion of the Uni- 
verse. Because of this, no matter when the change of co- 
ordinates is performed, when using FNC we cannot avoid 
a period when the evolution at second order in perturba- 
tions becomes relevant. To see this, considc;r an expand- 
ing universe with a single long-wavelength perturbation, 
in addition to short-wavelength modes as in Fig. 2. The 
metric in standard FNC takes a simple flat space form 
only on scales smaller than the Hubble scale or the 
physical wavelength of the long mode o/fci,, whichever 
is smaller. Therefore for the FNC metric to be valid in 
a region larger than the short modes, which is what we 
want to describe, we need ks aH (ks '> ki, being 
always true by assumption). But as soon as the short 
modes re-enter the horizon, ks ^ aH, they start evolv- 
ing, and acquire a non-trivial transfer function. In order 
to see that the evolution needs to be followed at non- 
linear order, consider our result for a general squeezed 
bispectrum Eq. (9). This expression contains terms in- 
volving the derivatives with respect to scale and time of 
the short mode power spectrum, which in general be- 
come non-trivial when ks ^ aH and the short modes 
start evolving with time. It is clear that these terms can- 
not be accounted for by just multiplying each of the per- 
turbations with the respective transfer function as linear 
perturbation theory would suggest. 

A nice way out of this complication is to use conformal 
Fermi Normal Coordinates (FNC), in which the metric 
is locally in the FLRW form, i.e. a homogeneous isotropic 
expanding Universe, rather than in the Minkowski form. 
These coordinates are valid up to a scale where fc^ 
denotes the wavenumber of the long-wavelength pertur- 
bation considered. At late times, the small-wavelength 
modes are well inside the horizon and have in general 
been subject to non-linear evolution. Then, transform- 
ing from the FNC frame to the usual FNC frame at a 
given spacetime point simply corresponds to a rescaling 
of the spatial coordinates. 

The final step is to relate quantities defined in the local 
FNC frame (e.g. temperature of the gas, or halo mass) to 
the photons actually measured on Earth. This involves 
the photon propagation ("projection") effects, such as 
lensing and redshift-space distortions [6, 11-14]. Here, 



we will adopt the "standard ruler" approach [6, 15] since 
it conveniently encompasses all projection effects for var- 
ious obscrvables. Specifically, the derivation in Schmidt 
and Jeong [6] assumes that a ruler, such as a ^(r) = const 
contour, corresponds to a fixed physical spatial scale on 
a constant-proper-timc hypcrsurface for a comoving ob- 
server, equivalent to a fixed spatial scale in FNC at emis- 
sion. 

Finally, at leading order, any physical (and necessarily 
non-gravitational) correlations between long-wavelength 
modes and small scale modes imprinted at early times 
and present in FNC then simply add to these projection 
effects. The framework discussed above, which we sum- 
marize in Fig. 1, thus connects the bispectrum calculated 
in any convenient gauge during inflation with late-time 
observables such as the CMB bispectrum or the large- 
scale clustering of tracers. It can be applied to any model, 
and we will give some examples for single clock infla- 
tion and the "fossil" scenarios studied in [16, 17]. We 
should mention that the importance of FNC in applica- 
tions to the calculation of inflationary perturbations has 
been recognized and crucially used previously in the lit- 
erature, e.g. in deriving various consistency relations [10], 
in proving that C, is constant at large distances at all loop 
orders [9] and in deriving the subleading corrections in 
ki/kg to Maldacena's consistency condition [8]. The issue 
of gauge artifacts in the squeezed limit of the bispectrum 
was also discussed in [18]. 

The outline of the paper is as follows: in Sec. II we de- 
rive some useful formulae to transform correlators from 
one set of coordinates to another. In Sec. Ill we in- 
troduce FNC, and show how the squeezed- limit three- 
point correlations calculated in the standard way trans- 
form to this frame. In Sec. IV we derive how correlations 
in FNC translate to observables, illustrating that, apart 
from projection effects such as lensing, the squeezed-limit 
three-point function in conformal Fermi coordinates is 
in fact the observed squeezed limit. In other words, if 
squeezed-limit correlations vanish in FNC, then the ob- 
served squeezed limit is only due to projection effects. We 
conclude in Sec. V and leave some technical details to the 
appendices. App. A and App. B derive the transforma- 
tion of the two and three point function in the squeezed 
limit, respectively, as one moves from one set of coordi- 
nates to another. In App. C we review the derivation of 
FNC and discuss their uniqueness in App. D. 
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II. TRANSFORMATION OF SMALL-SCALE 
CORRELATIONS AND THE SQUEEZED 
BISPECTRUM 

In this section we present some useful formulae to 
transform the two-point correlation function, Eq. (4), 
and three-point function and bispectrum in the squeezed 
limit, Eq. (9), from any one set of coordinates to another. 
The details of the derivation are left to appendix App. A 
and App. B. 

Consider a scalar field p{x), where x = (r, x) denotes 
the spacctime position. We work in comoving coordinates 
throughout, so that x° = t is the conformal time. Under 
a general coordinate transformation x — )• x'{x), the field 
transforms as 



p'{x') = p{x{x')) . 



(1) 



Consider a patch on a t' = const surface centered around 
the position Xq. We would like to derive the correlation 
^'(r',r') within that patch (defined with respect to the 
mean density /o'(t') over the patch) in terms of the two- 
point correlation function ^{r,T) of p in the unprimed 
coordinate system. Throughout the paper, we are consid- 
ering the case where the coordinate transformation varies 
slowly over the spatial patch, that is where we can ex- 
pand x in terms of x' as 



x\x') ^ xl, + A' ^ (^'■' -4') 
t{x')=to+Tj{x'^ -x2) , 

where x' = (r',x'), and 

xo =x(x'o,t') 
to=t{x'o,t')=t' + At 

dx'{x'o,T') 



(2) 



4' — 



dx'i 



(3) 



Note that a*^- and At are functions of Xq and r'. Here we 
have separated ^ into a zeroth and first order piece (the 
sign is chosen so that a'^ maps x^ to a;'* at linear order). 
The expansion to linear order in x' — Xg is sufficient, since 
we will only be interested in Fourier space contributions 
to the coordinate transformations up to order k\ (see 
App. A). As shown in App. A, the correlation function 
^'(r',r') of p' then becomes 

r (r', r') = [l - a,,- r'^c?^, + At 5,] ^(r', r') . (4) 

The dependence on x'^ of the time shift Eq. (2) does not 
contribute at leading order. 

We now turn to the squeezed limit of the three-point 
function {X' {x'^)S' {x[)S' {x'2)) , where S — Sp/p denotes 
the fractional perturbation to p (the results are identi- 
cal when considering 6p instead of 6). This limit corre- 
sponds to the case where |x3 — Xi| 3> |xi — X2I. Here, 
X stands for any other field {X does not have to be a 



scalar, although we suppress any tensor indices). Strictly 
speaking, X is to be understood as coarse-grained on a 
scale > |x2 — xi|. In the following, we will drop the 
primes on coordinates for clarity, since we will not em- 
ploy the unprimed coordinates anymore. Further, we 
adopt the notation X'(x) = X(x), that is, the long- 
wavelength (coarse-grained) field X is not modified un- 
der the small-scale coordinate transformation within the 
patch. In the squeezed limit, the three-point function 
describes the modulation of the local two-point function 
^'(|xi — X2I; Xq) at the location xq = (xi -|- X2)/2 by the 
long-wavelength field X evaluated at the distant point 
X3: 



squeezed 



(X'(X3)<5'(X1)5'(X2)) 

(X'(x3)e'(xi-X2;T)U'(xo)) 



(5) 



The precise location of the point xq in fact does not mat- 
ter in the squeezed limit (as we prove in App. B) , but we 
have chosen the midpoint as the most natural choice. 
We can now express {X'6'5') in terms of the three-point 
function in unprimed coordinates {X66) and contribu- 
tions from the transformation of small-scale correlations 
[Eq. (4)]. This yields 



(X'(X3)<5'(X1)5'(X2)) 

+ ^XAr(|x3 - Xo|)9-r 
+ (X(X3)5(X1)5(X2)) 



squeezed 



l-^Xa,,(|x3-Xo|)r'a^ 



^(xi -X2;t) 



(6) 



where ^xY{r) denotes the cross-correlation between X 
and Y = Uij, At. This expression gives the squeezed- 
limit three-point function of X' and 6' (in the primed 
coordinates) in terms of derivatives of ^(r, t) and the 
three-point function of X and d in the unprimed coor- 
dinate frame. As shown in App. B, we can derive an 
analogous expression in Fourier space, corresponding to 
the squeezed limit of the bispectrum: 



D n 1 1 \ squeezed 

5x'5'5'(kL,Ki,k2) = 



hi U 



kg d In ks 



P{ks,T) 



+ PxAr{kL)g^P{ks,T) 
5jf<5<5(kL,ki,k2) . 



ks=ki-|-ki/2 



(7) 



As before, Pxvik) denote cross-power spectra between 
X and Y, while P{k, t) denotes the power spectrum of 
p. We can now further decompose Q/j^ j clS 



riD 



(8) 



where n^, is the dimensionality of the space in which 
we define the correlations, and aj^ is traceless. Allowing 
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riD 7^ 3, in particular = 2 will become useful when 
dealing with projected observables such as the CMB. We 
then obtain 



1 1 \ squeezed 



Pxa{kL)P{ks,T) 



dln{k"s''P{ks,T)) 



din ks 



PxaT^kL) ,^2 



k^kidP{ks,T) 

kg d In ks 



+ PxArikL)-^Piks,T) 

-BxM(kL,ki,k2) . 



ks=ki+ki/2 

(9) 



Note that the contribution of the trace of aij scales as 
the logarithmic derivative of k'^^P{k), while the trace- 
free component afj couples to the logarithmic derivative 
of P(fc) itself. Instead of writing the bispectrum in terms 
of ks, one can also expand in fci/fci, yielding zeroth order 
terms (with kg — > fci) and first-order terms oc ki,/k\. The 
corresponding expression is given in App. B. 

The results of this section are useful when, say, Bxss is 
easy to calculate, while Bx'S'S' is more readily related to 
observations. This will be the case we encounter below. 



III. CONFORMAL FERMI NORMAL 
COORDINATES 

Our goal is to derive what the observable consequences 
arc (as opposed to coordinate artifacts) of a given primor- 
dial bispectrum. For this purpose we follow the steps in 
Fig. 1, starting with conformal Fermi Normal Coordi- 
nates (FNC). As we discussed in the introduction, if we 
want to avoid studying second order cosmological pertur- 
bation theory, it is not possible to directly connect the 
statistics of perturbations in the standard FNC frame of 
galaxies or the CMB (at the time when the photons we 
observe were emitted) with the statistics calculated when 
all perturbations are superhorizon. The reason is that the 
perturbative corrections to the Fermi frame metric are of 
order {HxpY, where xp are the physical FNC coordi- 
nates, so that the metric in FNC is not approximately 
flat on the scales of superhorizon perturbations. On the 
other hand, if we were to move from the global coordi- 
nates to FNC at a later time when the short wavelength 
modes have long re-entered the horizon, we would have 
to take into account their evolution after horizon entry 
at second order in the global coordinates. 

To overcome this obstacle, in this section we introduce 
a new coordinate frame which is valid throughout infla- 
tion and the hot big bang era (assuming nothing dramatic 
happens to superhorizon perturbations during this time) , 
and further connects easily with the FNC frame at late 
times. This provides a well-defined framework for study- 
ing observables related to primordial non-Gaussianities in 
the squeezed limit. To visualize the advantage of FNC, 



consider Fig. 2. Each one of the four panels shows the 
Hubble scale (solid circle) and the region of validity of 
the FNC (dashed line) at a different moment in time. 
The long (red line) and short (blue line) modes are also 
shown for comparison. One can immediately appreciate 
the fact that from very early time when the long mode is 
still inside (or about to leave) the horizon all the way un- 
til late time when the long mode has re-entered the hori- 
zon, the region of validity of FNC amply encompasses the 
short modes. This is in contrast with standard FNC for 
which the presence of the Hubble scale makes this im- 
possible. The primordial correlators, typically computed 
in comoving coordinates, take a simple form once both 
the short and long modes have left the horizon and cease 
to evolve. Any time after that moment and before the 
small-scale modes reenter the horizon constitutes a good 
time to transform from global conformal coordinates to 
FNC. 

Let us assume we are given the primordial correlations 
during inflation in some gauge, where the metric is writ- 
ten as 



ds^ = a^ir) [r?^^ -|- h^^] dx^dx" = a^(T)fif^^da;'*dx'' 



(10) 



Here we have denoted the conformal metric with a bar. 

In the absence of perturbations, h^^, = and g^^ = rj^i,. 
We will work up to linear order in metric perturbations. 

As in the usual Fermi normal coordinate construction, 
we consider a central timelike geodesic of a comoving ob- 
server. However, instead of constructing the FNC with 
respect to g^v, we construct conformal Fermi normal co- 
ordinates (FNC) with respect to g^^. That is, within a 
region around the central geodesic, the metric in these 
coordinates is approximately g^^ = a'^{TF)rjiiv The cor- 
rections to this metric, which determine the size of this 
region which we will call "FNC patch" , are given by sec- 
ond derivatives of the metric perturbations. Correspond- 
ingly, the size of the FNC patch is linked to the wave- 
length of the metric perturbations considered. In the 
following, we will be interested in (long-wavelength) per- 
turbations with wavenumber A;/,, so that the size of the 
FNC patch is of order k~[^. 

The transformation from coordinates to FNC, for 
a patch centered at the spatial origin at time x° = Tp is 
given by (see App. C 2) 



hoo{T)dT - (uj -I- hQi)x'' 



- -. [hoij + hoj,i + /ly x'x^ + 0[{x')^] (11) 



x%{x") =x''- v''{x^ - tf) + -h\x' 



x'x^ ^o\{xy\. 



(12) 
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FIG. 2: Illustration of the FNC patch throughout cosmic history, in comoving units. The wavy lines indicate perturbations. 
The solid circles denote the comoving horizon 1/aH, which coincides with the size of the usual Fermi coordinate patch during 
inflation. The dashed circles denote the size of the FNC patch, i.e. the region within which the metric is of the form g^^ — a^rj^i, 
with small corrections, a) During inflation, when the long wavelength perturbation /lij(kL) is generated (i.e. leaves the horizon), 
thin red line, b) Later on during inflation; hij (k_L) is far outside the horizon when the short wavelength modes J?(k) are generated 
(thick blue line). The perturbations are generated well within the FNC patch corresponding to the long- wavelength mode, c) 
Near the end of inflation. All perturbations are far outside the horizon. Nevertheless, the small-scale modes are still well within 
the FNC patch, d) At observation time, after the long-wavelength mode has reentered the horizon (solid circle indicating present 
horizon). The FNC patch now coincides with the usual Fermi coordinate patch, which is much larger than the small-scale modes 
which have been processed by nonlinear evolution since horizon entry (distorted thick wavy line). 



The inverse of this transformation is at linear order 
1 /"^^ 

x°{x'^) ^x°p + - / hoQ{T)dT + {v^ + hm)x'p 
^ Jo 

1 



x''{x"p)^V^{xl-TF 

1 r 

~ 4 



(13) 



x'px'p + 0[{x'pf] 



The metric in FNC then becomes 



-F 

9oo 



-F 
9^3 



1 

2 



dmdihoo + h'l^ 



x'px'P (14) 



9[rnKi] 



X pX p 



1 



-I- x'-px'^. 



Here, X^ab) = {Xab + ^ba)l'^, X\ah\ = {Xab - -^6a)/2, 

and primes denote derivatives with respect to t, while 
spatial derivatives are with respect to x*. The linear 
and quadratic coefficients in the coordinate transform 
Eq. (12) are uniquely determined by the requirement that 
the lowest order contribution to gf^^ — rjfii' is order (xp)^. 
However, given that the quadratic corrections are deter- 
mined by the cubic order terms in the coordinate trans- 
formation, one might wonder whether the quadratic cor- 
rections to g^^ are unique. One can show (see App. D) 
that gl^Q and are indeed unique. However, there is 
freedom to change the quadratic term in gfj if one allows 



the spatial coordinate lines to be non-geodesic at order 
{x'f. 

In summary, 



-F 

9^^ 



djh'^i,d,djh^^\x 



(15) 



Let us disregard the corrections cx /i'/^, djh'^^ for the mo- 
ment. The remaining corrections to the conformal Fermi 
frame metric scale as second derivatives of the metric per- 
turbation multiplied by the spatial coordinate squared. 
Thus, instead of being order {Hxp)^ as the usual FNC 
corrections , they are of order (fc^XF)^, where is 
the comoving wavenumber of the perturbation (in slight 
abuse of notation, we drop the bar over k since we will 
only be dealing with comoving Fourier wavenumbers) . 
Within the FNC patch of size ^ we can remove the 
effects of the long-mode perturbation (in comoving coor- 
dinates) at all times — through horizon exit and reentry. 

This very useful result only holds if h'-^, djh'^^ are of 
order didjhfj^u or smaller, which is the case for all models 
of single field inflation in which the background is an at- 
tractor solution. In principle one can engineer a model in 
which, for a short period of time, the background is not 
close to an attractor solution, but rather it is evolving 
towards one. In this case, the superhorizon modes can 
have, for a short period of time, sizable time derivatives 
[19, 20]. For example, if h'l^ is larger than didjh^^, cor- 
responding to significant superhorizon evolution of the 
spatial metric perturbation, then we cannot follow the 
conformal Fermi frame through the entire duration of in- 
flation. In case the evolving part of hij is isotropic (or 
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can be made so by a suitable change of coordinates), 

hij{0, t) = f{T)5ij + 0{didjh^^), (16) 

where denotes the origin in comoving coordinates 

around which wc construct the FNC frame, then we can 
absorb its effect into a modified scale factor, 



a(r) = a(T) 



1 - ifir) 



(17) 



In a modified FNC frame constructed with the scale fac- 
tor o(t) the rapidly evolving part now disappears. We 
thus see that in this case, a long-wavelength metric per- 
turbation does not simply shift the time and rescale spa- 
tial coordinates, but rather corresponds to a change in 
the entire background cosmology. This is possible only 
because the original "unperturbed" background cosmol- 
ogy, namely /(r) = 0, was unstable (i.e. not an attrac- 
tor). Given that the generation of small-scale pertur- 
bations depends on the background cosmology, we in 
general expect a non-trivial coupling between very long 
wavelength metric perturbations and small-scale modes 
in this case. If h"j is not isotropic, then we cannot absorb 
its effect into a modified scale factor. Instead, the differ- 
ent spatial coordinates are rescaled differently, leading to 
an anisotropically expanding Universe of Bianchi I type. 

Similar arguments apply to the case where djh'^^ is 
larger than didjh^v Here the lowest order effect will be 
a Bianchi I-type Universe, since isotropy is necessarily 
violated. 



With this, Eq. (4) yields 
C(f;f) = 



(20) 




haodr df 



e(f,t) 



where ^(f,f) denotes the correlation function of § in 
the local FNC frame. Now consider the bispectrum 
(X(kL)i?(ki)i?(k2)), where X is any perturbation. The 
results from Sec. II immediately show how this bispec- 
trum transforms into FNC. Denoting FNC quantities as 
X, 1^, we obtain 



squeezed 



-^PxH{k,)P.{ks) ^^^^^ 



+ \PxK 



kg d In ks 



+ Pxiir{kL)^P^{ks) 
^xtf^(kL,ki,k2;f) , 



ks=ki-|-ki,/2 

(21) 



where all correlations on the r.h.s. are evaluated at f . 
Here h = /i'j/3 is the trace of the spatial metric pertur- 
bation while hfj is the trace-free part. The significance of 
this result will become clear in the application to single- 
field inflation which we will consider next. 



A. Bispectra in global and FNC coordinates 

Let us assume that we can calculate the statistics of 
a scalar field i?(r, x) in the global coordinates Eq. (10). 
Specifically, we consider fluctuations at a given epoch 
during inflation, typically evaluated right after horizon 
crossing. For simplicity, we will restrict to a scalar field, 
although the extension to other spins is straightforward. 
Given the results from Sec. II and the coordinate trans- 
form Eqs. (12)-(13), we can then immediately derive the 
transformation of the two-point correlation ^.^(r, r) mea- 
sured in a patch around position Xq at conformal time r 
from the global coordinates to the FNC frame. Note that 
if ■& is non-Gaussian, then the correlation function mea- 
sured in a given patch will correlate with long- wavelength 
perturbations. 

In the present case, the primed coordinate system of 
Sec. II is Xp, while x'^ denote the global coordinates in 
the gauge chosen. We can then read off a^j and At from 
Eq. (13): 



At 



{T)dT . 



(18) 
(19) 



B. Single-fleld inflation in comoving gauge 

In the following, we restrict h^y to comoving gauge. In 
the notation of [4], we have 



/loo = - ZA^i 
hoi = Ni 

hij =2CSij + hJj{T,x), 



(22) 



where hX: is tranverse-traceless and contains the tensor 
perturbations. For the attractor solution of single-field 
inflation, the constraint equations in this gauge yield [4] 



H 



aH {any 



(23) 



To order k'j^/kg, we can thus neglect the contribution 
from the time shift At = — / NidT. Note that this is 
merely a consequence of the particular gauge chosen. The 
only remaining contribution to the transformation of the 
bispectrum [Eq. (21)] then comes from ( and the tensors 
h]]. For the scalar contribution, we obtain 



Bx^^(kLMM) = PxdkL)Pi){ks) 



3 + 



din P^ks) 
din ks 



-FBx,Jt?(kL,ki,k2), 



(24) 
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where throughout kg = ki + ki/2 = (ki — k2)/2. We 
are in particular interested in the bispectrum of the cur- 
vature perturbation i.e. iD = Since C is not a scalar, 
our derivation in Sec. II does not strictly apply. How- 
ever, it is straightforward to show that for the purposes 
of this transformation, the small-scale ( modes behave 
as a scalar. Recall that only the spatial transformation 
is relevant at the order we are interested. The comoving 
gauge condition, 5(1)^0 where S(j) is the inflaton pertur- 
bation is thus still satisfied in the FNC frame. Further, 
the spatial components of the metric transform as 



dx^ dx'' 



9ki{x)- 



(25) 



We now write Q{x) = Cl{x) + Cs{x), separating into 
long- and short-wavelength pieces on the scale of the 
patch within which the correlation function is measured. 
Then, the transformation to FNC removes Cl up to sec- 
ond derivatives, while (s is not affected since it does not 
contain any long-wavelength components. We thus ob- 
tain 



aXF)=Cs{x)+0{didjCL) 

Csixp) =Cs{x). 



(26) 



Thus, the short-wavelength perturbations transform ef- 
fectively as a scalar [Eq. (1)]. We obtain for the bispec- 
trum of curvature perturbations in the FNC frame: 



%-^-(ki,ki,k2) =Pc(fci)Pc(A;s) 



dln{klP^{ks)) 



dinks 

+ BfCc(kL,ki,k2), (27) 

where B^^^ is the three-point function calculated in co- 
moving gauge. As shown in [4, 21], this is in the squeezed 
limit given by 

Bccc(kL,ki,k2) = -[ns - l]Pc(A;i)Pc(fcs) , (28) 

where is defined through P({k) oc fc"^"*""'. Eq. (28) is 
usually referred to as the consistency relation. We now 
see that the first term in Eq. (27) exactly cancels the 
contribution from P^ff , leading to 



Pf-f^-(ki,ki,k2)=0(^ 



(29) 



We thus conclude that in single-field inflation, the bispec- 
trum in the squeezed limit is zero in the conformal Fermi 
frame, with corrections going as (kL/ks)"^- 

The case for non-scalar metric perturbations follows 
analogously. Decomposing the long-wavelength metric 
perturbation into polarization states. 



(30) 



we obtain 

+ BxM(kLMM)- (31) 

Assuming that the different polarization states are sta- 
tistically independent, we obtain for the tensor-scalar- 
scalar bispectrum [again using the transformation prop- 
erty Eq. (26)] 



1 
2 



lp,^ikMks)etj{^L)m^^^l§^ 
+ -B/(3fc(ki,ki,k2). 



(32) 



The squeezed-limit bispectrum in comoving gauge was 
also derived in [4]: 

P;,^Cc(kL,ki,k2) = -Ph,{kL)et^{kL)k'sk's^Pc{ks) 



1 d 



= -PhAkL)etj{]^L)k^sk's^J^/dks) ■ (33) 

We again see that the two terms in Eq. (32) cancel in 

single-field inflationary models. Thus, in single-field in- 
flation, the tensor-scalar-scalar bispectrum vanishes in 
the squeezed limit in FNC, so that there are no corre- 
lations between long-wavelength tensor m,odes and sm,a,ll- 
scale fluctuations in this frame. The lowest correction 
are again of order {kL/ks)"^- 

We can phrase the main result of this section as follows: 
at leading order, the squeezed-limit three-point correla- 
tions in single-field inflation, which obey the "consistency 
relation" , are equivalent to the statement that there is no 
correlation between infinitely long and short wavelength 
modes in the conformal Fermi frame, specifically, 



P;,.Kki,ki,k2) = 



(34) 



=+,x 



where h stands for any component of hij, and ks = |ki — 

k2|/2. 

Since this latter is a physical, gauge-invariant state- 
ment, it is expected to hold not only at leading order 
in spacetime perturbations but at higher orders as well. 
The FNC approach can also help elucidate why there is 
no such correlation up to order fc£ in single field models. 
One might wonder why such a correlation cannot be im- 
printed at early times when the long wavelength pertur- 
bation was inside the horizon, ki, ^ aH. Far inside the 
horizon, we can neglect gravity and are essentially dealing 
with a scalar field in vacuum which adiabatically tracks 
the slowly evolving background, and a given mode is only 
excited once its wavelength becomes of order the horizon. 
Contributions to correlations from within the horizon arc 
exponentially suppressed [9, 10] (see also [22]). On the 
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other hand, when ks ^ aH, the long-wavelength mode 
is far outside the horizon and its effects can be removed 
by a coordinate transformation up to order 

IV. CONNECTION TO LATE-TIME 
OBSERVATIONS 



As explained in Sec. Ill, we can follow the FNC patch 
all the way through the end of inflation and horizon re- 
entry of both long and short modes, provided that pertur- 
bations do not evolve significantly when they are outside 
the horizon. In this section, we show how the squeezed- 
limit correlations in the FNC frame can be related to 
observations made from Earth today. 

Let us assume we observe correlations of some field 
(5(x) at late times, i.e. at or after recombination. In 
linear theory, assuming adiabatic initial conditions, we 
can relate S to the curvature perturbation through some 
transfer function a, 

S{Kr) = a{k,T)C{k). (35) 

For matter density perturbations, a is defined in Eq. (71) 

below. Similarly, we assume the long-wavelength pertur- 
bation is evolved with some transfer function, 

hij{kL,T)=G{T)hij{kL), (36) 

where we have assumed the fci — >■ limit (the lowest or- 
der fci dependence of G(t) will be fc|; of course, G will 
differ for scalar and tensor perturbations) . As long as su- 
perhorizon perturbations evolve slowly (in the sense that 
h",djh' are smaller than or of order didjh), the confor- 
mal Fermi coordinate patch, i.e. the region over which 
corrections to the FNC metric are small, is essentially 
constant in (comoving) size throughout horizon exit and 
reentry of the short-wavelength modes. As we have seen, 
this applies in particular to single- field inflation models. 
Thus, at the conformal time Tem at which the photons 
we observe today were emitted, the metric in FNC is 

g^, = a''{T,^)r]^^ + 0{didjh) . (37) 

If we transform coordinates through 

X%= J\iT')dT' =t{T) 

x'p = a{T)x'p, (38) 
the metric in the coordinates {xp} becomes 

= V^.u + 0[{Hx'pf, (aaV)4]- (39) 

In other words, cLr6 the usual Fermi coordinates 

(FNC) defined around the same timclikc geodesic as the 
FNC. Depending on whether the long-wavelength modes 
for which we constructed the FNC patch have entered 
the horizon, either the order H or the order ddh^i, cor- 
rections will be dominant; however, this is not relevant 



for the discussion that follows. Eq. (38) corresponds to 
a rcscaling of the time coordinate, leaving the timelike 
unit vector and r = const hypcrsurfaces unchanged, and 
a time-dependent rescaling of the spatial coordinates. 
Since the spatial rescaling is the same everywhere on a 
x^p = const hypersurface, this implies that, for models 
that obey the consistency condition, Eq. (34) is still valid 
at late times in the FNC frame: 

Bh^s^5A^LMM\r,m) = o{^^ , (40) 

where hp, Sp denote perturbations in the FNC frame. 
Recall that the squeezed limit of three-point functions 
corresponds to the modulation of local two-point func- 
tions by long-wavelength perturbations. Thus, another 
way of phrasing this result is that a surface of constant 
correlation in the FNC frame, ^Sf{^f) = const, defines 
a standard ruler — a fixed spatial scale — as considered in 
[6], with corrections proportional to second derivatives 
of /lyii^ only. Hence, in standard single-field inflation and 
any other case where Eq. (34) holds, the ruler scale rp is 
statistically the same everywhere on a x'^p = const hyper- 
surface, and there is no correlation with long-wavelength 
perturbations. 

The apparent correlations induced between long- 
wavelength modes and small-scale correlations are then 
given by the ruler perturbations derived in [6]. We can 
use their results together with Sec. II to derive these con- 
tributions to observed squeezed-limit three-point func- 
tions, and make the connection with known results. 

In general, there are two effects modifying the observed 
two-point correlation within a given patch. First, there 
is the transformation from the local FNC Xp to the ob- 
served comoving coordinates x^, which arc inferred from 
the observed position n of the source in the sky and its 
redshift z through x° = tq — x(^)) = ^*x(^)) where 
X is the comoving distance-redshift relation in the back- 
groiuid (in case of the CMB, this is slightly modified, as 
we will discuss in Sec. IV A). Let us assume we observe 
a scalar field p{x) (whose perturbation is 6). Then, as 
described in Sec. II, p is given in terms of the field in the 
Fermi frame p as 

p{x) = p{xp{x)) . (41) 

At fixed observed redshift, we can write the transforma- 
tion from a;i? to i as 

aH 

= Ia'^x^p . (42) 

Here a=(l + z)^^,Tis the perturbation in proper time 
from a constant observed redshift surface [15], and A'j 
can be seen as the generalization to three dimensions 
of the magnification matrix. Note that T and A^j are 
gauge- invariant. Specifically, using the notation of [6], 
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Aij is given by 



(43) 



where V^^ = 5^^ — WW is the projection operator per- 
pendicular to the line of sight. C, Bi, and Aij are the 
gauge-invariant ruler perturbations derived in terms of 
the metric perturbations in [6]. For example, the trans- 
verse matrix Aij contains the magnification and shear. 

The second effect is a rescaling of the observed field by 
projection effects. If there are multiplicative projection 
effects, then we have 

/3obs(i) = [1 + c(£)] p{£) = [1 + c{i)] p{xf{5:)) . (44) 

Here we think of c as averaged over the patch within 
which we measure the correlation function of small-scale 
perturbations. For example, in case p is the number den- 
sity of some tracer, then projection effects such as grav- 
itational lensing modify the physical volume that corre- 
sponds to a fixed region in the observed coordinates x, 
thus rescaling the number density. Gravitational redshift 
and ISW effect lead to an analogous rescaling in case of 
the CMB. The factor 1 + c also rescales the fiuctuations 
in p within the region considered, and correspondingly 
rescales the correlation function by (1 + c)^, leading to 



{x)Fdi + -^T{i)dr + 2c{x) 
aH 



(45) 



where we have used Eq. (4) in the second line. Straight- 
forward application of the results of Sec. II then yields 
the bispectrum of the observed density perturbation 5 in 
the squeezed limit, in terms of the projection effects and 
the bispectrum of 6f in the FNC frame: 

-B|||(ki:,,ki,k2) = BspSpSf0^l,^i,^2) 

+ 2PspcikL)P5Aks) + -^P5MkL)^PsAks) 
aH OT 



+ 



idHklPsAks)) 



d hi ks 



P5Fa{kL)PiAks) 



+ Pspa-ikLykyki 



d In ks 



(46) 



where kg = ki -|- ki/2. The bispectrum of a two- 
dimensional projected field (in the fiat-sky limit) corre- 
spondingly becomes 

aH OT 



+ 



2 dlnis 



C5,a{£L)CsA^s) 



(47) 



where again £s = £1 + £l/2. 

The bispectrum BspSfSf in FNC frame is equivalent 
up to transfer functions (and radial projection, in the 2D 
case) to the bispectrum in FNC, which as we have seen 
vanishes in the squeezed limit for single-field infiation. 
In this case, the terms due to "projection effects" are the 
only remaining contributions. 

In the following we present two applications of 
Eqs. (46) (47): the CMB bispectrum in the squeezed 
limit, and the scale-dependent non-Gaussian halo bias. 



A. CMB bispectrum 

In this section, we illustrate the main result of the pre- 
vious section on the observed squeezed-limit bispectrum, 
Eq. (47), with the CMB. We will make direct connec- 
tion with the results of [13]. A more sophisticated and 
accurate treatment which is based on an closely related 
approach has been presented in [14]. 

We assume single-field infiation so that the FNC frame 
contribution vanishes. We will adopt the conformal- 
Newtonian gauge in this section. 



ds'^ = a^(r)[-(l + 2^)dT^ + (1 - 2*)dx2 



(48) 



although it is straightforward to derive the results in a 

general gauge. 

The observed CMB photons originate from the last 
scattering surface, which occurred at a fixed physical age 

of the Universe, that is, at constant proper time tp = 
t« for the comoving primordial plasma. The proper time 
of a comoving source passing through x at coordinate 
time T is given in the metric convention Eq. (48) by 



tF\r,^=i{T)+ [ $(x,T')a(r')rfT' 
Jo 



(49) 



where f(r) and f(t) are the physical-time - conformal 
time relations in the background. 

The value of is obtained by combining atomic 
physics with the mean observed temperature of the CMB 
today. The CMB temperature perturbations on scales 
that were super-horizon at recombination {£ < 100) orig- 
inate entirely from projection effects; in other words, the 
large-scale CMB temperature perturbations can be seen 
as a special case of the evolving ruler described in [15]. 
Essentially, the standard ruler is in this case given by 
the photon occupation number . In fact, since lens- 
ing conserves surface brightness, the only contribution 
to the fractional CMB temperature perturbation 0(n) is 
the redshift perturbation on a fixed proper time surface, 
which is equivalent to minus the quantity T defined in 
[15]: 



^ ' T d\na 



-r 



(50) 



where have used that T oi a ^ for a free-streaming black- 
body. As shown in [15], this immediately yields the CMB 
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temperature perturbation in the gauge Eq. (48) as 

e(n)=-iJ, / $[x,f(t')]dt' + iJo / <S?{0,f{t))dt 
Jo Jo 

rto 



1 

~ +HoJ $(0, f{t))dt - $o + vyo + 5i^s\ 

^0 



(51) 
(52) 



In the second Hne, we have assumed that the relevant 
perturbations are superhorizon, <I>(x) ss const (since we 
arc interested in ^ < 100 wlicirc thc^ ac;oustic contributions 
can be neglected) , and that recombination happened long 
after matter-radiation equality so that H^t^ = 2/3. 



1. Transforming from FNC to observer coordinates 

As we have seen, the CMB temperature on large scales 
is modified by the temperature perturbation Eq. (51), 
which is entirely due to the effects on photons as they 
propagate from the last scattering surface to the observer. 
Thus, the CMB temperature that would be measured in 
the local Fermi frame at emission is rescalcd by 1 + 0(n) 
in the observer frame, so that in the notation of Sec. IV, 
c(n) = e(n). _ 

The second ingredient is the coordinate transformation 
from FNC to the observer frame. Projected onto the sky, 
the matrix becomes the usual weak lensing distortion 
tensor Aij, which we can write to linear order as 



/ M/2 + J, 
\ 72 



72 

M/2- 



71 



(53) 



In particular, a = Ai. The magnification (a gauge- 
invariant quantity) is easily adapted from the results of 
[6] . The key difference is the meaning of the perturbation 
to the logarithm of the scale factor at emission, which 
here we call A In a,. By definition, 



Alna* = ^^-1. 



(54) 



Whereas the quantity A In o of [6] is derived for photons 
arriving with a fixed observed redshift, in the present 
case we need the corresponding expression for a constant 
proper time at emission. This is easily derived from the 
expression for the proper time Eq. (49). Requiring tp = 
and solving for r yields 



«(a;em) = a 



ft 



(l?(x, r) adr 



■■ a* X 



1-H, 



Jo 



X, r) adr 



(55) 



which leads to 
a* 



l-H, 



* r 

Jo 



X, T)adT , 



so that 



Alna, = -/f, / <b{x,T)adT . 
Jo 



(56) 



(57) 



We now use this result in the general-gauge expression 
Eq. (51) of [6], 



M 



2 A In a. 



1 



{h\ - /ill 



2k- 



-Axii 



(58) 



In Eq. (58), k is the coordinate convergence, i.e. the 
transverse (with respect to the line of sight) divergence 
of the transverse displacements, and Aa;|| is the displace- 
ment along the line of sight. In the metric Eq. (48) we 
have [6] 



dx— (x*-x)vi ($ + *) 



X* 



Aa;ii 



Alno* 



/ rfx[$ + *] T^Alna* 

Jo 

-I 



^{0,t)dt 



(59) 

(60) 
(61) 



where if, = i?(a*), and we have made the same assump- 
tions about superhorizon perturbations and matter dom- 
ination as in Eq. (51). The magnification then becomes 



^ = 3* 



(-24-) 



2k Aa;|| 

X* 



10^ 2 ^ 

= — $ + 2k Axu 

3 X* 



(62) 



Here we have used that ^' = $ during matter domination. 

Using the relation between convergence k and shear 7 
in ^-space, we then have in the flat-sky approximation 

CTc{e)=CTT{i) 

in 2 

Craii) = CTMit) = ^CT4.e„ W + 2Ct«(^) - — CtAx,, 

•-> X* 

(63) 
(64) 



In Eq. (63), we have approximated C'T<i>em (■^) /3 as the 
CMB temperature power spectrum without ISW contri- 
bution, and neglected the contribution from the linc-of- 
sight displacement Axy . The latter is small due to cance- 
lation along the line of sight except for the very smallest 
£ for the time delay [first term in Eq. (60)], and sup- 
pression by (if»x*)~^ ~ Vl089 in case of the second 
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term of Eq. (60) (the third term in Eq. (60) only con- 
tributes to the monopole). See [14] for a quantitative 
cvahiation of these contributions. Finally, for the contri- 
bution from the time shift, Ctt{^l) = —CTri^L) since 
Q = —T- However, the CMB power spectrum observed 
today Ctt{£s) only evolves on the Hubble scale today, 
so that 



^CTT{iL)^CTT{is) 

aH ^ ' OT 



^CTT{iL)CTT{is) . (65) 

aH 



Thus, this contribution is suppressed by ^ l/\/l089 with 
respect to the leading contributions, and we will neglect 
it in what follows. 



2. Squeezed-limit CMB bispectrum 
Inserting these results into Eq. (47) then yields 

B{£l,£i,£2) = 2CTc{h)CTT{es) 

1 a In 4Ctt(^s) 



+ CTa{eL)CTTi£s) 



2 dlnis 



+ COS20 CTii{iL)'^^^^^^^^\cTT{is), (66) 



where £3 = ^z,+-^i/2, and cos^ = t-L-i-s- Apart from the 
very small contribution from Axy , this is an exact result 

for {\\c squecizcid-limit CMB bispectrum as long as the 
approximation of the Sachs- Wolfe limit is accurate for 
modes of wavenumber {(.l ^ 100). Note that no such 
assumption has been made about the short-wavelength 
modes £s- The cross-correlation between CMB temper- 
ature and k can further be decomposed as 

CTfi{£L) = Ic^^^ai^L) + Cs.,,^i,{£L), (67) 



i.e. an early-time correlation with the potential "I>cm) and 
a late-time correlation between k and the ISW contribu- 
tion. The latter contribution has recently been detected 
by the Planck satellite [23]. 

In order to compare with Creminelli et al. [13], we 
neglect the ISW-lensing correlation, and use the result of 
[12] for early-time correlations in the Sachs- Wolfe regime. 



(68) 



= 2Ctt{£l) - CMiL) 
- 6 cos 261 CttC^l) 



9 In 4Ctt(^s) 



d\n£s 
d\nCTT{(-s) 



ain^s 

Ctt(^l)Ctt(^s) (1 + 6 cos 20) 



Ctt{(-s) 
, d\nt%CTTm 



d\n£s 



(69) 



Further, we neglect the distinction between Ctt{£l) and 



This agrees with the final result of [13]. 

One can show (App. C of [6]) that the contribution 
of metric perturbations with wavenumber k to the CMB 
temperature as well as the roller perturbations Ad , 7 scale 
as (k/Ho)'^ in the limit A; — > 0. The contributions in the 
low-fc limit can be interpreted as the lowest order correc- 
tions to our local conformal Fermi patch. Unless there 
is some physical coupling between currently superhori- 
zon and subhorizon modes [24], the observable imprint 
of any superhorizon perturbation (scalar or tensor) is 
suppressed by {k/Ho)^. 



B. Non- Gaussian halo bias 

For Gaussian initial conditions, the distribution of 
large-scale structure tracers (such as galaxies, clusters, 
etc.) in the large-scale limit follows the distribution of 
matter. More accurately, this holds on a constant proper 
time slice [7, 25] . While the abundance of tracers depends 
on the amplitude and shape of local small-scale fluctua- 
tions, in the Gaussian case these are statistically the same 
everywhere. That is, at fixed proper time in their respec- 
tive FNC frame, all observers sec statistically the same 
small-scale fluctuations. These then do not contribute to 
correlations in the tracer abundance on large scales. Non- 
Gaussianity in the primordial perturbations can however 
couple small-scale fluctuations to large-scale perturba- 
tions. As we discussed in Sec. II, this in fact precisely 
corresponds to the squeezed limit of (at lowest order) 
the three-point function. Thus, neglecting for the time 
being any projection effects in going from the FNC frame 
to the observed positions and redshifts of tracers, there 
is a non-Gaussian scale-dependent bias if and only if the 
amplitude of small-scale fluctuations in the FNC frame 
correlates with long-wavelength perturbations — that is, 
if there is a non-zero squeezed-limit bispectrum BspSpSp' 
where 5p denotes the matter density perturbation in the 
FNC frame. 

Using Eqs. (35)-(36), the bispectrum oi Sp is re- 
lated to the bispectrum of curvature perturbations in the 
FNC frame through 

)a(fc2,Tem) 
X %-^-(ki,ki,k2), (70) 
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where 



a{k) = 



5 2 k^T{k)g{z) 
3 3 0„if2(l + ^) 



(71) 



is the relation in Fourier space between the density and 
the curvature perturbation, T{k) is the matter transfer 
function normalized to unity as fe — > 0, and g{z) is the 
linear growth rate of the gravitational potential normal- 
ized to unity during the matter dominated epoch. 

In order to derive the scale-dependent bias, we now as- 
sume that the tracer abundance is mostly sensitive to the 
variance of the density field on a scale R, typically chosen 
to correspond to the Lagrangian scale of the tracer. See 
[25, 26] for a more general and detailed discussion. The 
contribution to the tracer two-point function is then pro- 
portional to ((5i?(0)(5F,_R)^(r)), where 6f,r denotes the 
FNC-frame density field smoothed on scale R. As shown 
in, e.g. [27, 28], the large-scale scale-dependent bias is 
then given by 



Abiki^) = ^^a-\kL)AaUkL) 



/J3 L 
^a(A=i)a(|k + k: 



B.-.-.-(A:L,A:i,|k-Fki| 



XWR{kl)WR{\\^ + \^ 



(72) 



where Wji{k) is the filter function in Fourier space, and 
we have dropped the Tcm arguments for brevity. As long 
as kL is much less than the typical values of ki in the 
integrand (fci 0.01 /i Mpc~ ), the bispectrum here is 
evaluated in the squeezed limit [29]. Let us perform a 
leading order expansion in this limit, 

%-f(fci,fci, |k + kil) = (^/o + /2^) P^{kL)P^iks), 

where ks is as defined before. Here, we have assumed 
a scale-invariant bispectrum (otherwise, the coefficients 
/o: /2 are in general functions of ks), and no dependence 
on the angle between k^ and ks as is typically the case 
(but see [30] for a counterexample). The lowest order 
piece oc /o leads to a scale dependence of 



Ab{kL) oc fo a-\kL) oc fo kl^T-\kL) 



(73) 



which corresponds to the usual scale-dependent bias from 
local non-Gaussianity. The second part leads to a scale- 
dependence of 



Ah{kL) (X h kl a-\kL) « /2 T-\kL) 



(74) 



which is much weaker and only relevant on scales fci > 
0.02/1 Mpc~^ where the transfer function departs from 
unity. Since fo = in single field inflation [Eq. (34)], 
we conclude there is no scale-dependent bias induced on 
scales kL < 0.02 /i Mpc~^ in these models. 

The scale-dependent bias discussed here refers to the 
FNC frame of the tracers; there are contributions to the 



observed scale-dependent bias from projection effects, i.e. 
from transforming from the tracer FNC frame to ob- 
served positions and redshifts (in our own FNC frame). 
Those have been derived in [7, 31-34], but are generally 
quite small. Matching to the scale-dependent bias from 
local non-Gaussianity, they correspond to /nl < 1 for a 
wide range of tracer parameters [34]. 



V. CONCLUSIONS 

In this work we have presented a simple and complete 

framework to translate the squeezed limit of any pri- 
mordial three-point function into observables such as the 
squeezed limit of the CMB temperature bispectrum and 
the scale-dependent halo bias. In Fig. 1 we have spelled 
out the various steps of the computation and the respec- 
tive choices of gauge. The different gauges have been 
chosen in such a way that the whole computation from 
horizon exit during inflation until observation on earth 
can be performed using only linear perturbation theory. 
This required the introduction of a conformal version of 
the well known Fermi Normal Coordinates (used in [9, 10] 
for inflationary correlators), in which the spacetime is lo- 
cally FLRW as opposed to Minkowski. As an example 
we have applied our formalism to standard single-field 
slow-roll inflation. We have shown that Maldacena's con- 
sistency condition [4] in the squeezed limit of the scalar 
bispectrum implies that the signal in the CMB bispec- 
trum and halo bias vanishes exactly. Although this re- 
sult was already known, our approach provides a simple, 
concise and physically clear derivation. This calculation 
can straightforwardly be generalized to higher N-point 
functions in the squeezed limit. 

In addition, our approach shows that, for exactly the 
same reason, the consistency relation for the scalar- 
scalar-tensor bispectrum implies the absence of any 
observable correlation between long-wavelength tensor 
modes and short-wavelength scalar perturbations ("fos- 
sils") as proposed for example in [16] (see also [35, 36]). 
Taking the expression for the tensor-scalar-scalar bispec- 
trum in comoving gauge Eq. (33) at face value, Ma- 
sui and Pen [16] showed that the signal-to-noise ratio 
squared for detecting a scale-invariant spectrum of ten- 
sor modes with wavenumbers in the range [fcrnimfemax] 
scales as k^^^ (if k^ain is horizon scale or smaller). On 
the other hand, we have shown that in single-field infla- 
tion the only contributions to the observed correlations 
between long-wavelength tensor modes and small-scale 
perturbations are due to projection effects, which scale as 
{k/aHY times the tensor perturbation where aH is to be 
evaluated at the time of emission [37, 38]. In this case, the 
signal-to- noise ratio squared roughly scales as fcmax) and 
is thus dominated by the smallest observable wavelengths 
rather than the largest. Note that due to the oscillatory 
nature of the tensor mode transfer function (i.e. the fact 
that tensor modes propagate), the observable effects of 
tensor modes with k ^ aH are usually highly suppressed 
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[38, 39]. On the other hand, just as for the scalar bis- 
pectmm, a detection of a squeezed-hmit tensor-scalar- 
scalar bispectrum through the measurements described 
in [16, 17], at a level larger than expected from projec- 
tion effects, would rule out single-field inflation. 



The framework we have developed in this paper is also 

useful in other classes of models that have not yet been 
directly related to observations. One example is resonant 
non-Gaussianity [40-44] , which is a generic prediction of 
many models of axion inflation (see [45] and references 
therein), in particular inflation from axion monodromy 
[46-48]. In these models the inflationary potential has 
sinusoidal modulations, which lead to oscillations in the 
time evolution of the background. Since these oscillations 
average to zero over a period, they can give corrections to 
the slow-roll parameter that are larger than is usually al- 
lowed and nevertheless be perfectly compatible with cur- 
rent power spectrum data (sec e.g. [47, 49, 50]). The reso- 
nant bispectrum satisfies Maldacona's consistency condi- 
tion [42, 51], but because of oscillations the amplitude 
of the primordial bispectrum in the squeezed limit is 
not suppressed by small slow-roll corrections. In light 
of this, one might wonder whether these models lead to 
some detectable signal in the scale dependent halo bias. 
Cyr-Racine and Schmidt [52] showed that indeed reso- 
nant non-Gaussianity produces oscillations in the mass 
dependence of the non-Gaussian halo bias, which is a very 
unique signature. However, the leading contribution in 
the large-scale limit to the effect derived in that paper is 
only a coordinate artefact as shown here. While the halo 
bias thus has to asymptote to a scale-independent value 
in the large-scale limit, we expect some interesting effects 
on intermediate scales. To compute the actual size of this 
effect it is important to re-write the primordial resonant 
bispectrum in terms of FNC. It would be very interesting 
to perform this analysis using the framework constructed 
in this paper. There arc also other models of the early 
universe (both infiationary and not) in which the metric 
perturbations do not freeze outside of the horizon. For 
example, it has been argued in [19, 20] (sec also [53]) that 
this allows one to violate Maldacena's consistency condi- 
tion. Another interesting model with a peculiar behavior 
in the squeezed limit is Khronon infiation, studied in [54]. 
It would be interesting to use our approach to derive the 
observational prediction of these and similar models. 



On the other hand, multifield inflationary models in 
general feature a non-trivial bispectrum in the squeezed 
limit. In this case, curvature perturbations evolve outside 
the horizon, and the small-scale fluctuations are sensitive 
to their presence as they essentially evolve in a different 
FRW background (or, more generally, in a homogeneous 
anisotropic Universe). This reasoning can be used to con- 
nect the squeezed-limit correlators in multifield inflation 
to late-time observables in a similar way as outlined here 
for single-field inflation. 
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Appendix A: Coordinate transformation of 
two-point correlations 

Consider a scalar field p{x) = p(x, r). Under a general 
coordinate transformation x — x'{x), p transforms as 



p'{x') = p{x{x')) . 



(Al) 



We define the unequal-time two-point correlation of p, 
measured within a patch p through 

(p(x,ri)p(x + r,T2))p 

= ( Wl) + ^^(X, Tl)] [p{T2) + 6p{x + r, T2)] )p 

= p(Tl)p(T2)+e(r,Ti,T2), (A2) 

where a subscript p indicates an average over the patch, 
and we have defined p(t) = {p{x,T))p, so that by con- 
struction the average over the patch of 6p{x, t) vanishes. 
The correlation function introduced above is then defined 



as 



^(r,Ti,r2) = ((5p(x,Ti)5p(x-|-r,T2))p. 



(A3) 



We will consider an alternative definition, the correlation 
function of 5p/ p, below and show that it transforms in 
the same way. Further, we will denote the equal time cor- 
relation function as ^(r, r) in the following for simplicity. 

Now let us consider the equal-time two-point correla- 
tion for p' , measured in the same patch: 

(p'(x;, r')/>'(x'2, r'))p = (p(xi, ri)p(x2, r^))^ , (A4) 

where X2 — x'^ = r', and we have defined for conve- 
nience Xi = x(x'i,r'), X2 = x(x2,t'), and analogously 
for Tl , r2 . By assumption, the coordinate transform is 
slowly varying over the patch. We then Taylor expand 
the coordinate transformations around the center point 
x(, = (x;+x'2)/2. 



Tl = To 



1 5x'(x'o,tO . 
2 dx'i 

2 dxi 



"0 



dx'3 



T2 



°^2 dxi 



(A5) 
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where Xq = x(xq, t'), and tq = t(xq, t'). Wc will discuss 
below why going to linear order in r' is sufficient. We 
then have 

(p'(x;,r')p'(x'2,r'))p 

, iar(x'o,r') a iar(x;„r') d 



dT2 



Tl=T2=To 

(A6) 



since ^(r, ri, is symmetric in ri and T2. Thus, at linear 
order in r', the equal time correlator in primed coordi- 
nates is an equal time correlator in unprimed coordinates. 
We thus have 



(p'(x', r')p'(x' + r', r'))p = p'(ro) + ^^v' , ro) 



4' — 



dx'3 



(A7) 



We now define the correlation function of p' in the same 
way as that for p [Eq. (A2)], which yields 

(p'(x', t')p'(x' + r', r'))p = p'^t') + ^'(r', r') , (A8) 

where 

p'(t') = (p'(x', t'))p = (p(a;(x', t')))p = p(to) . (A9) 

The second equality again holds to linear order in x' — Xg, 
at which order the average of r over the patch at fixed r' 
is just Tq. We will discuss this definition of ~p' in App. A 1 

below. Thus, the first term in each of Eqs. (A7)-(A8) 
agrees, and hence so must the second term: 

^'(r',r')=e(x2-xi,ro)=^(Ar',T(x'o,T')) • (AlO) 

We can similarly derive what happens to ^ = ■f/p^, which 
is the correlation function defined in terms of 5p/~p. In 
this case, we have 



Thus, we see that ^ transforms in the same way as ^, 

namely 

|'(r',T') = |(X2 -Xi,To) - |(Ar',T(x^,r')) . (A14) 

There are thus two ways in which the general affine 
coordinate transform affects the correlation function in 
primed coordinates: first, there is the spatial transfor- 
mation of the separation vector, r — ^ Ay'; second, the 
primed correlation function is evaluated at a different 
point in space and time (here, wc have only made the 
time shift explicit in the notation, but note that the cor- 
rrelation function on the right hand side is to be evalu- 
ated within a patch centered on x(xq, t')). 

The linear order expansion in Eq. (A5) neglects higher 
derivative terms in the coordinate transformation, i.e. it 
is valid in the limit that the transformation is slowly 
varying over the patch. The terms we are neglecting 
in Eq. (A5) correspond to terms of order fc| and higher, 
where fci, is the wavenumber of the long- wavelength mode 
that contributes to the coordinate transformation. Since 
we are neglecting terms of order k\ throughout this pa- 
per, it is sufficient to expand to linear order. 

A second simplification occurs because we are only con- 
sidering three-point correlations in this paper. In this 
case, it is sufficient to consider the linear response of 
^'(r', r') to the coordinate transformation. We thus write 



AT = r(x^„T')-r' 



(A15) 



and work to linear order in ■ and At. Further, we can 

neglect the effect of the spatial shift in the position at 
which ^ is evaluated. This is because the contribution 
from this shift, 

e'(r',r') D [x(x^,r')-x;,] . V,,,e(r',r')U^ (A16) 

is second order, since both the spatial shift and the lo- 
cation dependence of ^ are at least linear order in the 
long-wavelength perturbations. Eq. (AlO) then simpli- 
fies to Eq. (4), 

C'(r'; t') =\l- aijr"dl, + At C(r'; r') . (A17) 



(p(x, r)p(x + r, r)) = f{T) [l + |(r, r) 
For p', this yields 

(p'(x',t')p'(x' + r',T'))p =p'(to) [1 +|(x2 -Xi,To) 

(A12) 

Now we define the correlation function of p' in the same 
way as that for p [Eq. (A2)], which yields 

(p'(x',r')p'(x' + r',r'))p =p'^(r') \l + i'{r',T') 



(All) 

1. Coordinate transformations in the presence of 
non-trivial backgrounds 

There is a conceptual point to clarify about the trans- 
formation properties of perturbations under changes of 
coordinates. There are two different but equivalent ways 
to work with coordinate transformations when study- 
ing perturbations on top of some non-trivial background. 
The first one consists in changing every field in the same 
way, whether it is going to be treated as a background 
or not in the rest of the computation. This "democratic" 

approach is conceptually very simple. For example, if a 
(A13) 

scalar field transforms as p'{x') = p{x{x')), so does its 



15 



background p'{x') = {p'{x')) = {p{x)) = p{x). Using 
this point of view, perturbation of the scalar field around 
that same background must transform in the same way 
as well. To see this, consider 5{x) = p{x) — p{x) (notice 
that {5) = 0). Then 

5'{x') = p'{x') - -p'{t') = p{x) Pit) = S{x) , (A18) 

In the second approach to the problem one splits the 
field in perturbation and background in such a way that 
the latter is invariant under coordinate transformations, 
i.e. p'{x) = p{x) (in the cosmological context, the back- 
ground quantity is only a function of t) . This simplifica- 
tion comes at the cost of more complicated transforma- 
tion laws for the perturbations, which now do not trans- 
form as a scalar field. Considering the same example as 
above and a transformation x' = x + e, at linear order in 
e one gets 

5'{x') = p'{x') p'{x') = p'{x') p{x') (A19) 
= p{x) - p{x) + e^d^,p{x) = 5(x) + e^9^p(x) . 

Notice that in both approaches p'{x') = p{x{x')), while 
the two different conventions above concern only the 



transformations of background and perturbations. A 
similar discussion can be given for tensor fields such as 
the metric. Again one has two choices: either one defines 
metric perturbations that transform like the component 
of a tensor, in which case the metric background changes 
when one changes coordinates or the background is as- 
sumed fixed in any coordinates and metric perturbations 
have additional terms in their transformations. 



We warn the reader that in the derivation of App. A 

we have used the "democratic convention" as opposed 
to the fixed background convention which is standard in 
the cosmological literature (e.g. [55]). However, it should 
be kept in mind that aroimd a homogeneous background 
with a long wavelength perturbation, the additional con- 
tribution in the last term of Eq. (A19) is just a con- 
stant or a pure gradient at the order we are working in, 
and hence does not contribute to the correlation function 
evaluated on much smaller scales r. This means that our 
final result Eq. (A17) is valid independently of the con- 
vention used for defining the perturbations. 



Appendix B: Squeezed-limit three-point function from transformed two-point correlation 

This section derives the squeezed-limit three-point function and bispectrum from a coordinate transformation of 
the two-point function given by Eq. (A17), 

^'(r; r) = [1 - a^.r^S^ + At 9,] C(r; r) , (Bl) 

where we have dropped the prime on coordinates since we will only deal with primed coordinates in this section. 
Specifically, we want to derive the three point function (X'(x3)5'(xi)(5'(x2)) in the limit where |x3 — Xi| |xi — X2I 
(squeezed limit). Here, X and 6 stand for any perturbation variables with mean zero (of course, the auto-three-point 
function X = 5 is a special case). In this limit, the three-point function quantifies the modulation of the local two-point 
function ^'(|xi — X2I) by a long-wavelength perturbation X' (modes with wavelength much less than |x3 — Xi| will 
not contribute to this correlation). As discussed in Sec. II, the coordinate transformation only acts on the small-scale 
fluctuations, so that we set -^'(x) = X{x). Thus, 

(X'(x3)<5'(xi)(5'(x2)) ^''"^^^^X'(x3)e'(xi -X2;r)U,(,^)), (B2) 

where we have defined the point of evaluation of the long-wavelength X perturbation as 

Xc = cxi -I- (1 - c)x2 . (B3) 

Choosing x^ along the axis connecting Xi and X2 is necessary since Eq. (B2) describes a homogeneous and isotropic 
three-point function. The midpoint Xq = (xi -I- X2)/2 corresponds to c = 1/2. We will see that the choice of c does 
not influence the flnal result below. Using Eq. (Bl), this yields 

(X'(X3)^'(X1)5'(X2)) [-Cxa.,(|x3 - X^Dr'^^ + ^xAr(|x3 - ^c\) Or] ^(xi -X2;r) + (X(x3)^(xi - X2;T)U(xe)) 

= Hxa,, (|X3 - ^clVdi + ^XAr (|X3 - X^j) dr] ^(xi - X2; t) + (X(x3)^(xi)5(x2)) . (B4) 

This expression gives the squeezed-limit three-point function of X' and S' (in the primed coordinates) in terms of 

derivatives of ^(r, r) and the three-point function of X and 5 in the unprimed coordinate frame. 

The left-hand side of Eq. (B4) is clearly symmetric under Xi 0x2. For a general choice of Xc, this does not hold 
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for the right-hand side, so we should symmetrize: 



(X'(x3)5'(xi)y(x2)) 



squcczci 



d 1 



+ ^XAr(|x3 - Xc|) dr + ^XAr(|x3 - x'^|) dj 



^(xi-X2;r) + (X(x3)<5(xi)^(x2)) 



(B5) 



where x^ = (1 — c)xi + cx2. Of course, if we choose x^ to be the midpoint of xi and X2, the two permutations are 
identical. 

We now derive the Fourier-space analog of Eq. (B5). In terms of the cross-power spectra -Pxay (^)) PxArik) and 
the auto power spectrum Ps{k), we have 



(X(x3)a,,-(x,)) = I g|P;,„^^.(fc^)e^'^-[-3-xc] 
(X(x3)Ar(x,)) = J |^PxA.(A;i)e^''-[— =1 



(B6) 



The bispectrum in the squeezed limit fci, <^ ki, k2 becomes 

(X'(ki)5'(ki)5'(k2)) = j rf3x3 j d^xi y d3x2e-*(''-^3+''^^^+''^^^)(X(x3)<5'(xi)(5'(x2)) 



squeezed 1 

" 2' 



Pxa^A^L)^ [khPsiks)] + PxAr{kL)^P5{ks) 



+ {c^l-c) 



- ks=ki-|-cki. 



X i2nfSD (ki + k2 + ki) + (X(ki)(5(ki)5(k2)) . 



(B7) 



Thus, the transformed bispectrum has the proper delta function ensuring the triangle condition, and we can identify 
the 5'— bispectrum in the squeezed Hmit as 



-Bx'5'<5'(kL,ki,k2) 



squeezed 1 | p^^^^^j^^^^^ [k^Ps{ks)] + PxA. (^l) (^s) 

-BxM(kL,ki,k2) 



+ (c ^ 1 - c) 

ks=ki+cki, > 



PxaJkL) U, + j Ps{ks) + PxAAkL)-^P6iks) 



+ Bxss{^L,'k-i,'k-2) ■ 
Since we are working in the squeezed hmit, we can expand in ki/kL- We have 

F(|ki + ck^l) = P(fci) + cki • k^l^i -^F(fcO + 0{kl/kl) 



+ (c ^ 1 - c) ^ 

ks=ki+cki, ) 

(B8) 



ki dlnfci ' 



4-2 i-2 



kl 



dPsik.) _dPsik,) +ek,.k,^^^P,(fc0 + O(fciAf). 



d In fcg d In fci 



fci (i(lnA;i)2 



(B9) 
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We already see that all corrections arc linear in c, so that the two permutations c, 1 — c add up to 1. Thus, the result 
becomes independent of the choice of c, and the squeezed limit bispectrum in terms is, to order {k\/k\)^ 



kf d In ki 



din All 



—rr-i ^ Pxat{kl) 

dmki 



dr dlnfci j 
(BIO) 



Since k2 = — ki +0{kL/ki), this expression can be equivalently written in terms of k2 instead of ki. Since the choice 
of c is arbitrary, we will use the most natural choice, c = 1/2, in which case Eq. (B8) becomes simply 



_ , , , squeezed 

-Bx'5'5'(kl,,ki,k2) = 



PxaJkL) i^Sij + j Psiks) + PxAr{kL)^Ps{ks 

+ Sx55(kL,ki,k2) . 



ks=ki+ki,/2 



(Bll) 



While this expression is equivalent to Eq. (BIO) up to order fc£, we will work with this result since it is more compact 
and convenient. 

Including the effect of a mean density modulated by X following Eq. (44) in Sec. IV is now an obvious generalization. 
We obtain 



squeezed 



Bx'S'5' (ki,,ki,k2) 

+ Bjc55(ki,ki,k2). 
We can now further decompose Uij as 



Pxa,,{kL) Uj + j Pgi^ks) + Px Ar {kL)-^Ps{ks) + 2Pxc{kL)Ps{ks) 



J ks=ki+ki,/2 

(B12) 



no 



a6ij +aip 



(B13) 



where Hd is the dimensionality of the space in which we define the correlations, and afj is traceless. In particular 
no = 3 for three-dimensional observables such as galaxy densities or 21cm flux, and n/j = 2 for projected quantities 
on the sky such as the CMB. We then obtain 



-Bx'5'5'(kL,ki,k2) 



squeezed 



1 d 
ud dinks 



PxaikL) (1 + ^ ) Psiks) + PxaL (^l) ,2 



fc^4 dPsjks) 

kg d In ks 



d 



+ PxAr{kL)^Ps{ks) + 2Pxc{kL)Ps{ks) 

Pxa{kL)Ps{ks)-'^ 



d\n{k^s- Ps{ks)) 



d\n kc 



+ PxaTXkL) 



+ Bx55(kL,ki,k2) 

ks=ki+ki,/2 

44 dPs{ks) 



kg d In ks 



d 



+ PxAr{kL)Q^Ps(ks) + 2Pxc{kL)Ps(ks) 



Bx55{^lM,^2) . (B14) 



ks=ki+ki,/2 



Note that the contribution of the trace of Ojj scales as the logarithmic deerivative of k'^°Ps{k), while the trace- free 
component aj, couples to the logarithmic derivative of Ps{k) itself. 

I 



Appendix C: Fermi Normal Coordinates 
1. FNC construction and metric 



do lead to a metric of the form 



In this section, we show how for a general metric, the 
coordinates defined through Eq. (C2) below in fact 



(Cl) 
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To quadratic order, the Fermi coordinate transformation 

is given by 



1, 



x'^{xp) ^ — x'^ {P) + {ei)p Xp 2^0:^ p^^i')p^^:)'^P^F'''F 

(C2) 



tp 



Straightforward algebra yields the transformation matrix 
as 

dx° 



dx'^ 



ie^rp + A^x'p + Oix'p), 



(C3) 



where a subscript P denotes the evaluation at the point 
on the central geodesic specified by Xp. Here, 



A^^ = ieofp[dpie,r] 



- /37 



(C4) 



Using the fact that the unit vectors (e^)" are assumed 
to be parallel-transported along the central geodesic, we 
have 

which can be used to bring A'^j into a uniform expression: 



A" 



■01 



(C5) 



Next, we expand 

gcp{x{xF)) = gap 
Finally, using that 



+ {e3)p[9'ygafi]pXp. (C6) 



(C7) 



by construction of the orthonormal tetrad (note that this 
holds for any velocity Vi at order v), we obtain 

9„Axf) = gc.p{x{XF))Q^Q^ 

F F 

+ 0{xl). (C8) 
Assuming that we have a metric connection so that 

V^5„^ = d^gap - 2r^(„5A/3) = 0, (C9) 

we see that the order xfp correction to the metric van- 
ishes. 

So far, we have not used the condition that the central 
curve whose tangent vector is (eo)^ is a geodesic. How- 
ever, if we want the metric in FNC to be of the form 
Eq. (CI) at more than just one point along the central 
curve, we clearly need at lowest order 

(eo)W^ [ga0{e^r{e.f] p = 0. (CIO) 
This in particular implies 

{eoV^^ie^r = (Cll) 

along the curve, which is precisely the condition for the 
central curve to be a geodesic. 



2. Transformation from general coordinates to FNC 

We now explicitly derive the transformation Eqs. (12)- 
(13) into the FNC frame. We write the conformal metric 

as 



(C12) 



and work to linear order in /i^y. The tetrad around P is 

given by 



(eo)" = ( 1 + ^/loo, V 



V, + /in,-, ^/ 



(C13) 



Note that we treat %; as 0{h) here, as we arc assuming 
comoving observers. Since {ca)^ = '^,^+C>{h) and F^^ = 

0{h), we have for the transformation into FNC at lowest 
order in h: 



x'^{x"p)=P'^ + {e,rpX^p-^-r^^ 



VpX'^p , 



(C14) 



where P^ = x^{P) are the global coordinates of the cen- 
tral geodesic, so that = dP^ /dx^ and we define analo- 
gously P^ = Xp{P) for FNC. The conformal proper time 
P° defines the time coordinate of FNC, namely P° = Xp. 
Without loss of generality, we can choose the spatial ori- 
gin so that P'^irp) = at some fixed proper time Tp. 
To lowest order in h and v, the conformal proper time is 
related to the global time by 

X°{P)=X%+1 [ " /loo(0, T)dT . (C15) 

^ Jo 



The Christoffel symbols are given by 

rr,- = I + h% - K,'-] , (CI6) 

where indices are raised and lowered with r/^^. Thus, 

1 



1 



/ioo(r)(ir + {vi + hoi)xp 



4 ^^-^ 



1 



pjf'p 



'■Xf+ ^ I hoo{T)dT + {vi + hoi)xp 
^ Jo 



1 



[hoij + hoj^i + h'^j] XpX^p 



x\x%)=P''--h\x\ 



(C17) 

ipX'^p ■ 
(C18) 



Here all terms linear in perturbations, namely h^i, and 
v^, should be evaluated along the central geodesic P^. 

For x'^ Tp, the spatial position of the geodesic differs 
from the origin in the x'^ coordinate system by an amount 
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Ax'^v'', which is Hnear in perturbations (we will shortly Here. S^^ij in general depends on the affinc parameter 



find from the geodesic equation). Therefore, up to 
terms quadratic in metric perturbations, we can evaluate 
h^^, and at the origin of the FNC x'p = P'{x%J = 0. 

At linear order in h, the inverse transformation is very 
easily derived: 



1 

2 70 

[hoi,j 
1 



hoQ{T)dT - {vi + hoi)x^ 



x%{x°') = x" + -h%x' + 



1 
4 



X-X-' 



(C19) 



,fc 



x'x-' 



(C20) 



Note that when hoi 7^ there is no solution in which 

the observer is at rest. We now focus on this case now, 
where we can set hij = Hqq = 0, since their effect is 
simply additive and was considered above. For hio 7^ 0, 
one needs to know the relation between and hoi that is 
enforced by the geodesic equation for P. Assuming as we 
did previously that the velocity = dP^/dx^ is of the 
same order as hoi, the geodesic equations at linear order 
in perturbations in global coordinates are given by 



d 



P" = 0, 



\dx° ) 



fdP°Y _ 

00 Va^,' 



where we used the fact that T^^ is at least linear in met- 
ric perturbations. The first equation tells us that proper 
time coincides with the global time coordinate at this or- 
der, up to two arbitrary integration constants, which we 
fix by choosing P° = x^. This is no surprise since we are 
assuming /iqo = 0. Then, from the second equation, we 
find v'' = dP' /dx'^ — —hoi, where we have set another in- 
tegration constant to zero following the assumption that 
doP"^ is linear in hoi- Specifically, for Xp close to tf, we 
have 



x\P) = v\sPp 



TF 



(C21) 



so that to leading order 



x^{x%) = v^{x% - tf) +x%- ^h%x'p 

1 

~ 4 

x^ix") =x^ - v''{x° - tf) + \h\x 



i ^3 

pXp 



(C22) 



x'x-' . 



along the central geodesic, i.e. the Fermi-frame time coor- 
dinate tp- Let us consider a general inertial frame {.t'^} 
constructed around point P. Without loss of generality, 
we let P be at the origin of both the {x'^} and {xp} 
coordinate systems. The requirement that {a;'^} be in- 
ertial, i.e. that g'^^^ = 1}^^ at P and that d'^g'^^^ = at 
P, restricts the relation between the coordinates to be at 
least of cubic order: 

x% = x'" - \c"0^sx'^x'^x'^ + 0{x'*). (D2) 



Here, 



which implies 



C" 



\J tJU pi JU pi tLi pi 



(D3) 



(D4) 



i.e. the last equality states that C'^p^g — ^ sj^- 
metric in the primed frame then becomes to order x'^: 

9',xv = V1J.U + [Sij.u-y6 - 2C(^^)^5] x''^x'^. (D5) 
By assumption, S^^^s obeys the "FNC condition" 

Sfii^os = = Sf^iyjo- (D6) 

Imposing the same restriction on g'^,^ leads to the condi- 
tion 

C(Aii^)0i5 = = C(^^)^o- (D7) 
Using the symmetry properties Eq. (D4), we obtain 



(m0)7<5 



0. 



(D8) 



This says that there is no freedom in the order ©(xl^) cor- 
rections to the Fermi-frame metric components gQQ and 
for coordinates satisfying the condition Eq. (Dl) (but 
see below). However, there is some freedom in choosing 
the spatial part gfj, since a non-zero Cijki is allowed by 
Eq. (Dl). 

Finally, there is an additional freedom in the choice of 
coordinates which does not affect the metric at order x%. 



(C23) Namely, we can choose C^Q^^ij ^ 0. Specifically, 



(D9) 



Appendix D: Uniqueness of Fermi Normal 
Coordinates 



In this section we investigate what residual coordinate 
freedom remains when requiring that the metric be of 

"FNC form" : 



F 



, F'^ F' 



(Dl) 



where we have introduced a 3-vector ak and a sym- 



metric 3-tensor Aij. By construction, C( 



(0m)7<5 



0. 



This corresponds to a coordinate transform of (note that 



CO 



kij 



-Oik Aij) 



X'° = X 



x'^ = Xp — a'^x 



F CtkX F A.ijXpX'p 

Xp AijX pX^p. 



(DIO) 



20 



In other words, this coordinate transform is a spatially 
location-dependent Lorentz boost with a velocity given 

by 

v^{x^p) = a'' AijXpX^p, (Dll) 



which vanishes quadratically on the central geodesic. 
This type of coordinate transform leaves the Fermi frame 
metric entirely invariant at order x'^. 
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